Scalar and vector Keldysh models in the time domain 
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The exactly solvable Keldysh model of disordered electron system in a random scattering field 
with extremely long correlation length is converted to the time-dependent model with extremely 
long relaxation. The dynamical problem is solved for the ensemble of two-level systems (TLS) with 
fluctuating well depths having the discrete Zi symmetry. It is shown also that the symmetric TLS 
with fluctuating barrier transparency may be described in terms of the planar Keldysh model with 
dime-dependent random planar rotations in xy plane having continuous 5*0(2) symmetry. The case 
of simultaneous fluctuations of the well depth and barrier transparency is subject to non-abelian 
algebra. Application of this model to description of dynamic fluctuations in quantum dots and 
optical lattices is discussed. 
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The model with infinite correlation range of fluctuat- 
ing fields V(r) proposed by L.V. Keldyshi is one of few 
exactly solvable problems in the theory of disordered elec- 
tron systems. The approximation 



D(r - r') = (V(r)V(r')) = W 2 



(1) 



makes identical all diagrams for the electron Green func- 
tion (GF) in the order V 2n . As a result summation of 
diagrammatic series in the " cross technique"— reduces to 
the problem of calculation of combinatoric coefficient A n 
(number of pairwise coupling of scattering vertices). In 
fact An = (2n — 1)!! is the total number of identical 
diagrams in the order 2n. The perturbation series is 
summed exactly^, and one deals with averaging of en- 
semble of samples with constant V but the magnitude 
of this field randomly changes from realization to real- 
ization. In momentum space the correlation function ((T|) 
transforms into -D(q) = (2n) 3 W 2 S(c[). The electron GF 
in Keldysh model is averaged with Gaussian distribution 
function characterized by the variance W 2 . 

This model is not widely used in current literature 
because it is difficult to propose an experimental de- 
vice, where the conjecture ([T]) could be realized (sec, 
however^). In the present paper we discuss a realiza- 
tion of Keldysh model in time domain. In this case the 
analog of infinite spatial correlation is the long mem- 
ory effect, which can be realized in many physical situa- 
tions (see below). The structure of perturbation series in 
time-dependent Keldysh model (TDKM) is the same as 
in original one, and the long characteristic times of dy- 
namical correlation play the same part as infinite range 
spatial correlation of static random potentials. Since the 
time axis is the only coordinate in this problem, its effec- 
tive dimension is "0+1". Moreover, the TDKM admits 
natural generalization of original Keldysh model. We will 
show that the dynamical fluctuations in time domain may 
be both of scalar and of vector character. The kinematic 
constraint existing in the vector TDKM results in elim- 
ination of essential part of diagrams in cross technique, 
but the summation of perturbation series is still exact. It 



results in 2D Gaussian averaging for the GF in dynamical 
random field. 

Leaving for the last section the discussion of real sys- 
tems, where TDKM arises as a description of generic 
disorder, we start with a toy model of an ensemble of 
non-interacting two-level systems (TLS) in a randomly 
fluctuating environment. In standard realization of TLS, 
namely a double- valley well, particles arc distinguished 
not only by conventional quantum numbers but also by 
their position in the well characterized by the index 
j = l,r of the left (I) or right (r) valley. The barrier 
between the valleys is characterized by the tunneling ma- 
trix element Ao. The Hamiltonian of isolated TLS has 
the form 



H. 



(0) 
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Utf)-A (4c r + K.c). (2) 



Here rij = CjCj is the particle occupation number, Sj is 
the discrete energy level in the valley j and U is the in- 
teraction parameter for two particles in the same valley. 
The condition U 3> A is usually assumed. We consider 
spinlcss particles, having in mind that the theory can be 
applied both to interacting bosons and fermions (elec- 
trons) with frozen spin degrees of freedom. To be specific 
we discuss tunneling electrons as an example. 

We start with the singly occupied TLS, where the 
constraint N = J^i Tii = 1 is imposed on the Hamilto- 
nian, introduce pseudospin operators cr + = cjc r ,a~ = 
n r , and reduce ^ to 



4c;,2cr 2 = m 



H- 



(0) 



-S a z - A a x - fx (N - 1). 



(3) 



in the pseudospin subspace. Here the asymmetry param- 
eter Sq = e r — £1 play the role of effective "magnetic" 
field, the Lagrange parameter fj,Q controls constraint. 

The scalar fluctuation field is introduced as random 
fluctuations of TLS asymmetry, namely as a time depen- 
dent field 6p(t)=So+h p (t) determined by its moments 



h p (t) = 0, hp{t)h p {t + r) = D{r). 



(4) 
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Here the ovcrlinc stands for the ensemble average. Thus 
we reduced the original model to the effective spin Hamil- 
tonian in magnetic field with random time-dependent 
component. The problem can be reformulated as a 
study of propagation of fermions along the time axis in 
the presence of time-dependent random scalar potential 
5 p (t)(n r — n/)/2, and the cross technique may be used in 
calculation of the propagators^. 

The analog of Kcldysh conjecture in time domain is 
a slowly varying random field ^exp(— 7t). A very long 
relaxation time T re ;~ 1/7 with small 7 is presumed, so 
that the noise correlation function is given by 



D(tu) = lim 

7^0 LU 



2C 2 7 



2tt( 2 6{lu) 



(5) 



(the noise correlation function (|5|) is normalized in such 
a way that the corresponding vertices are dimcnsionless). 
In this limit the averaged spin propagator describes the 
ensemble of states with a field i5 = const in a given 
state, but this constant is random in each realization. 
The " planar- type" TDKM may be derived in a similar 
way. For this sake one should introduce a random com- 
ponent in the tunneling matrix element, A = Ao + A p (t) 
with A p (t) = and make similar conjecture (JSJ) about 
the correlation function F(t) = A p (t) A p (t + r), namely 
approximate its Fourier transform by 



F{u) 



lim 

7^0 LU 



4£ 2 7 
" +7 2 



^ 2 5(tu) 



(0) 



Thus, we treat our toy Hamiltonian ([3]) in the following 
way. First we study the limiting case So » Ao, were the 
interdot tunneling is considered as a perturbation to the 
longitudinal term affected by stochastization in the scalar 
TDKM. This problem may be solved by means of the 
standard techniquei^ generalized for the time-domain 
case. Then we turn to the case of transversal random 
tunneling potential and solve this problem by means of 
correspondingly modified "planar" TDKM. 

We start with the scalar TDKM and treat the term 
if 11°' = So(n r — ni)/2 as a zero order approximation with 
electron occupying the level e\ in the ground state. With- 
out stochastic perturbation the role of the tunneling term 
Hj_ = — Aqo~ x is in admixing a charge transfer exciton 
to the ground state with the corresponding energy level 
shifts, Ei^ r — > e^ r =F Aq/5 for the ground and excited 
states, respectively. The time-dependent perturbation 
stochastisizes this simple picture. 

Let us introduce the retarded propagators for the 
scalar TDKM 

Gf, s (t f ) = (c J (t)c](t')) R = -i{[ Cj (t)c](t')} + ). (7) 

and consider their evolution on the time axis under the 
influence of random component hj(t) ("random longitu- 
dinal magnetic field" in pseudospin notation), first as- 
suming 7— >0, A0/7 — > 0. After averaging G^(t — t 1 ) in 



accordance with (|4|) and making the Fourier transforma- 
tion by means of ([5|), we come to the series 



GfAe) = 9j {e) 



(8) 



Here gi(e) = (e + in)^ 1 and g r = (e + 5q + iri)^ 1 are the 
bare propagators, A n = (2n — 1)!! is the above mentioned 
combinatoric coefficient (see Fig. [IJ where several first 
irreducible diagrams are shown). 





FIG. 1: Irreducible Feynman diagrams for scalar TDQM. 
Solid and wavy lines stand for gi(e) and D(u). 

Like in the real space Keldysh models, the series ([5]) 
may be summed by means of the integral representation^ 
for (2n— 1)!! . Then changing the order of summation and 
integration (Borel summation), one comes to the follow- 
ing equation for the left valley GF 



1 



2 /2C 2 



dz 



VI) 



(9) 



Remarkably, the single electron GF in this model has 
no poles, singularities or branch cuts. Similar procedure 
may be applied to the Green function G^. As a result 
of this Gaussian averaging the " magnetization" a z is re- 
duced and the corresponding response to transversal field 
is modified accordingly. 

Next, we formulate the "planar" TDKM for the bare 
Hamiltonian H$ s ((2) with impenetrable barrier Ao — > 
and symmetric valleys, ei = e r , isolated from each other. 
The transverse random perturbation is introduced by 



H p {t) = A p (t)a+ + A*(t)a 



(10) 



so that the inter-valley tunneling is stochastisized by 
means of averaging in time-domain with correlation func- 
tion (JB]) . The tunnel matrix clement A is transformed as 



A, 



-<pi) 



under the gauge transformation Cj 



Cje lVj , and we presume A p to be a complex variable. 

Unlike the scalar TDKM, the noticeable part of dia- 
grams in the perturbation series disappears due to the 
kinematical restrictions a + a + = o~ a~ =0 (see also^). 
Only the diagrams with pseudospin operators ordered as 
. . . o~ + o~ cr + o~ . . . survive in the expansion for the GF of 
planar model 

00 

Gf: p (e) = 9j (e) + B n (V20 2n gf l+1 (e). (11) 

n=l 

The vertices in the cross technique are now "colored" in 
accordance with two terms entering the random poten- 
tial. The vertices with different colors have to be ordered 
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FIG. 2: First non- vanishing vertex correction to the Green 
function self-energies in the planar Keldysh model. Black 
and white sites correspond to two terms in the Hamiltonian 
(??); transversal pseudospin correlation functions F(iv) © 
are represented by dashed lines. 



in alternating way, and the correlation lines connect only 
the vertices of opposite color (see Fig. [2J) 

As a result of the above kinematic restrictions, the 
combinatoric coefficient B n = n\. Then we use the inte- 
gral representation^ for n\ and transform the series (fTTj) 
into 



Gf >p {e)= gi {e) 1 + 2^ 



tdt 



(Figs. [U [2J and corresponding skeleton diagrams^ are 
still valid. However, the important reservation should 
be made: the self energy cannot be treated as a renor- 
malization of bare pole because the bare and dressed GF 
are connected by non-local integral operators [see Eqs. 
([5]).(|13p]. Nevertheless, the ordinary differential equa- 
tion connecting the GF and its derivative over energy 
can be obtained for both versions of TDKM. This equa- 
tion was found for the scalar Keldysh model in». Here 
we derive this equation without appealing to the Ward 
identity and then generalize the derivation procedure for 
the planar model. 

To calculate the derivative dG(e)/de for the scalar 
TDKM, we start with expansion (index s is omit- 
ted below for the sake of brevity), ft is convenient to 
count the energy off the position of chemical potential 
Ho [see Eq. ([3])] in the middle between the levels e^ r . So, 
we shift the energies e — * e± = e =F So/2 for the left and 
right GR, respectively. 

The same procedure, which leads to Eq. © for GR 
gives the following equation for its derivative: 



Here we substituted t 2 for the variable z. Then changing 
the order of summation and integration, we transform 
G^(e) into the integral 



ffj( e ) 



G j>P (e)-J Q 2tdt - 1 _ 2t 2 eg 2 {e y 



(12) 



Taking into account the explicit form of the free propaga- 
tor gj(e), we change the integration variable once more, 
t = u/^/2^, an d transform (|T2"]) into 



Jo 



50 udu 

w 



1 



1 



e — u + it] e + u + irj 



(13) 



Now we introduce the "cartesian" coordinates, x = 
u cos 0, y = usin0, so that u = \/ x 2 + y 2 and dxdy = 
ududtfi. The angle independent integral (fTB")) may be 
rewritten as 

nR _ 1 f+°° dxe-^l^ f+°° dye-y*' 2 ? 
Gj,p{£) - 2 



1 



\/x 2 + y 2 +irj e+ \j x 2 + y 2 + if] 



(14) 



This result is a natural generalization of the one- 
dimensional Gaussian averaging ((9|) characteristic for 
the scalar TDKM to the two-dimensional Gaussian av- 
eraging of planar random field with purely transversal 
{xy) fluctuations. Only the modulus of random field 
r = \J x 2 + y 2 is averaged, whereas the angular variable 
remains irrelevant due to the in-plane isotropy of the sys- 
tem. Like in the scalar model, the averaged GF has no 
singularities. 

Although the GF lost its pole structure, the stan- 
dard Feynman rules for construction of irreducible parts 



dG 

de 



3 2 (£a) 



CV2TT 



2C 2 



C 2 £c 



IT] 



dz. 



(a 



Calculating the integral and substituting 



g (Sa) = £ a we come to the differential equation 



c a dg(sq) 
de 



1 - e a G{e a ) 



(15) 



similar to that obtained for the real space scalar Keldysh 
model^i. 

Generalization of this procedure for the symmetric pla- 
nar model (5q = 0) is more cumbersome. We start with 
differentiating the series (|11[) over the energy. Then the 
analog of Eq. (fT2")) for the derivative has the form 



dG 
de 



-9 2 {e) 



1 



2tdt 



2(2t 2 - l)t 2 £V(e) 
1 -2t 2 £ 2 g 2 (e) 



The subsequent variable change which gave Eq. (|13|) for 
the GF gives for its derivative the following equation 



dG 

de 



-9 2 {e) 



1 



1 fJi J2 



2 u 4 e 



(16) 



where 



Jn = J™ dzz n eiq ? (-7^) W - z ) + + z)] 

After some manipulations, these integrals are represented 
via the GF for the planar model (fl3|) : 



J 2 = 2e( 2 G -2^ 2 , J 4 = -4£ 4 + e 2 J 2 



(17) 



Substituting these integrals in Eq. (|16p . we come even- 
tually to the differential equation 



dG{e) 
de 



1 - eG(e) ( 1 - § 



(18) 
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which is obviously the generalization of Eq. (|15|) . These 
two equation may be rewritten in a unified way: 



&aG a ^ s 
EGr, 



1 = eo 2 p 



(19) 



The solutions of these equations satisfying the boundary 
condition G(e — > oo) = e _1 is given by ([9l [T3|) . It is worth 
noting that the differential operator in the right hand side 
of the second equation is nothing but div e in polar coor- 
dinates. This form reflects effective two-dimensionality of 
Gaussian averaging in the planar TDKM. Now one may 
introduce vertex parts using the analogy with the Ward 
identities 



de 



g: 



ld_ 

e de 



(eG- 1 ). 



(20) 



These vertices, together with equations (|19[) will be use- 
ful for calculation of response functions of our TLS (see 
below). Like in the self energy parts (Figs. [HO, the pla- 
nar TDKM lacks most of diagrams of scalar model due 
to the kinematic restrictions: the sites in the vertices of 
triangle are of the same color, black and white sites al- 
ternate, and dashed lines connect sites of opposite colors. 
First nonvanishing vertices for both models are shown in 
Fig. El 
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FIG. 4: Density of states in scalar (left) and planar (right) 
model (all units are arbitrary). 



double hump Gaussian structure. The information about 
position of electron in right or left valley is completely 
lost at C/^o > 1- hi the planar model for symmetric 
TLS instead two-peak structure due to avoided crossing 
characteristic for coherent tunneling, we get a pseudogap 
around zero energy due to stochastic tunneling, which 
survives at any variance £. 

Let us allow now the fluctuations of both longitudinal 
and transverse components which corresponds to simul- 
taneous fluctuations of the well depth and barrier trans- 
parency for the case of symmetric TLS (Sq = 0). We 
solve this problem by means of path integral formalism. 
The Lagrangian and corresponding action are defined on 
the Keldysh contour K (see e.gJ^-) 



C{t) 



cAdtCj — H, S 



K 



£{t)dt. (23) 



A" 



Here Cj , Cj arc Grassmann variables describing the elec- 
tron. The time-dependent gauge transformation Cj(t) — > 
Cj(t)e lipj ^ , converts the fluctuation of the well depth to 
the fluctuation of the phase of the tunnel matrix clement 
under the choice 



FIG. 3: First non- vanishing vertex diagrams for the F in scalar 
(left) and planar (right) TDKM. 

The density of states (DoS) in stochastisized TLS is 
given by the imaginary parts of GF ©, (O for scalar 
and planar TDKM, respectively: 



v s {e) 



cxp 



1 



cosh Ke 



2|e|exp 



2£ s 



(21) 
(22) 



In the scalar model iy s (e) is a superposition of two Gaus- 
sians centered around Ei and e r , respectively. In the pla- 
nar model Vp(e) is represented by a single Gaussian with 
a dip "burnt" around zero energy (Fig. [4j. 

Switching on the tunneling term H± in scalar model 
without random field, we come to the picture of two lev- 
els mutually repulsed due to coherent interdot tunnel- 
ing and broadened due to incoherent time-dependent in- 
tradot fluctuations. The spectrum is still gapful at small 
enough ratio (/Sq, If the C < <5o, the DoS merges into 



hj(t')dt'. 



(24) 



We therefore identify the longitudinal and transverse 
noise with phase fluctuations of the barrier transparency 
and fluctuations of the modulus of the tunnel matrix el- 
ement, respectively and unify them in the path integral 
description. 

The ensemble averaging 

(...)„oiae= / dh p P[{h p ) J dA* p dA p P tr (A* p ,A p )... 

is done with the help of probability distribution functions 
for longitudinal and transverse fluctuations 



Pi = 



^= exp ( - / dtdt'hp^D- 1 ^ - t')h p (t')) 
\/2n \ Jk J 



Pf.r — 



2< 



■ exp 



K 



dtdt'A* p {t)F- 1 {t-t')Ap{t') 



The GF's can be calculated by means of generating func- 
tional corresponding to Keldysh action Sk hi a standard 



wa y 12 ) 13 . In the "infinite memory" limit ([5]), ([6]) we easily 
express the GF of the electron in the symmetric double 
well potential 



1 



C£ 2 (2tt) 3 / 2 
dw*dwe-^ 2/2 ^ 



O 

e±z 



(e + it]) 2 — z 2 — | -lc | 2 



(25) 



Noticing that the GFs do not depend on the well index 
j and performing the integration over angles in spherical 
coordinate system we get 



erf 



1 

2i 



dpp exp 



1 



2£ 2 
1 



p + ir\ e + p + irj 



(26) 



The Eqs fl23 ESI) generalize © and dT4j) for the 
case of anisotropic vector Keldysh model. The three- 
dimensional Gaussian averaging in (|25p stands for vector 
character of the random field distributed on an ellipsoid. 
Typical size of semi-axes is defined by the variances of 
longitudinal and transverse noises. Like in the scalar an 
planar models, the averaged GF has no singularities. The 
angle (f> dependence is absent in (|25[) due to the in-plane 
isotropy of the model Pt r = Pt r (|A| 2 ) which is preserved 
here. The limits of strong easy axis £— >0 and easy plane 
£^0 anisotropy correspond to scalar and planar models, 
correspondingly. 

We notice that the DoS for vector TDKM is also 
characterized by a pscudogap. The energy dependence 
v v (e)~£ 2 for small energies should be contrasted with 
v s (e)~const for scalar v p {e)^e for planar model behav- 
ior. This dependence reflects the fact that the probability 
to remain close to initial level position decreases with in- 
crease of effective dimensionality of the problem due to 
stochastization of the complex tunneling matrix element. 

One should specially mention the degenerate case of 
the isotropic vector TDKM, characterized by £=£=A and 
describing rotation of pscudospin on a Bloch sphere. Per- 
forming simple algebra we get an equation for GF 



eG v - 1 



A G„, 



and the Ward Identity 



r "-?& (£2 ^ 1) ' 



(27) 



(28) 



The DoS for the isotropic model is given by the expression 



v v (e) 



1 



A 3 V2V 



e exp 



2A 2 J 



(29) 



Our next task is calculation of response func- 
tions. In scalar TDKM the longitudinal suscepti- 
bility is given by the correlation function X||( w ) = 





FIG. 5: Diagrams for bare loop (transverse susceptibility) and 
loop with dressed vertex (longitudinal susceptibility). 



i J dt exp (iujt)(a z (t)a z (0)}ii, which is represented by two 
loops with vertex corrections (Fig. [5]). Here both solid 
lines correspond either to j = I or to j = r. We confine 
ourselves with calculation of static susceptibility, u> — > 0. 
In order to work at finite temperatures we turn to Mat- 
subara Green functions functions Q(ie n ) (similar calcula- 
tion can be done on the Keldysh contour in the real time 
path integral formalism) and susceptibility 

)Qj(ie n )T(ie n ,iuj m ) (30) 

Then the first of Ward identities ([20)) provides us with 
the exact equation for the vertex 



C 2 (G*) 2 r« (e,0) = eG&-l 



(31) 



giving an access to the exact evaluation of x(0)- Com- 
bining ([3H)l with (|31[) and analytical continuation of ©, 
we find 



xn(o) = -£ 



ydye-y 2 / 2 f (2y-aS )C \ 

~v^r nF { 2T ) (32) 



Here np(x) is the Fermi distribution function. The 
asymptotic behavior of static susceptibility x(0) is 



X(0) 



1/T, T»(C,5 ) 
1/C, C»(T,<5 ) 



(33) 



There is no vertex correction to transverse susceptibility 
X-l(0) = (er+,cr-)_R. It is given by the bare loop formed 
by the left and right GF. In case of symmetric TQD 
(do = 0) and big variance C 3> A this function is as 
smooth as X||(0) with changing £ and T. Its asymptotic 
behavior is given by the same Eq. (|33p . The physical 
sense of these results in obvious: the information about 
position of electron in a given well is lost at strong enough 
stochastization £ ^ T. 

Next, we calculate the susceptibility for planar model 
in case of symmetric TLS with 5q = 0. In this case 
X_l(0) differs from X||(0) by factor 2, so that it is enough 
to calculate the latter one. Now we appeal to the second 
equation from ([19]) with T p defined in f|20[) . Then like 
in previous case the calculation of X||(0) is reduced to 
finding the combination ie n Qp{i^n) — 1- Straightforward 
computation gives 



t Jo 



y 2 dye-y 2 ^t & nhf^ 



(34) 



The behavior of X|| as a function of T and £ is close to that 
for scalar model, including the asymptotic dependence 
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(|33|) . The equations for static susceptibilities in isotropic 
vector TDKM can be easily found with help of ([25 1 125 ]) 

^) = ^/°°^ 2/2 tanh(g). (35) 

These susceptibilities are characterized by the same 
asymptotic behavior as those for scalar/planar models. 

The toy model of noninteracting TLS under dynami- 
cal stochastization demonstrate some generic properties 
of TDKM: (i) the loss of characteristic spin or pseudospin 
behavior at variance exceeding temperature; (ii) the ef- 
fective two-dimensionality of Gaussian averaging in pla- 
nar TDKM as its main distinction from scalar model; 
(iii) the effective three-dimensionality of vector TDKM. 
These features survive also in more realistic situations. 
One of possible applications of this theory is the problem 
of electron tunneling through double quantum dot in a 
regime of strong Coulomb blockade, where the source of 
stochastization is a random time-dependent gate voltage 
applied to one of the valleys^. The case of N = 2 was 
considered, where the starting Hamiltonian H^°> is that 
of Eq. (J2j with added spin index. In this case the scalar 
TDKM may be used in the limit of slow fluctuations ([5j, 
the double quantum dot looses its spin characteristics at 
low T, and the Kondo-type zero bias anomaly is smeared 
accordingly^. The important difference between spinful 
and spinlcss TLS models is in their symmetry. The sym- 
metry of TLS with N=2 considered in£ is 50(5), which is 
reduced to 50(3) for low-energy part of excitation spec- 
trum, so the Lie algebra is non-abclian. However, it was 
possible to introduce scalar TDKM due to abelian char- 
acter of time-dependent random gauge field. 

The symmetry of the scalar TLS with N=l is given by 
the discrete group Z 2 with abelian algebra. The gauge 
transformation allows to identify the fluctuations of the 
scalar model as U(l) fluctuations of the phase of tunnel 



matrix clement. In the planar model one deals only with 
the planar (xy) rotations, so the relevant symmetry is 
50(2) with still abelian algebra. The phase fluctuation in 
that case are discrete Z% (0=0, 7r to provide the condition 
A p =0), while modulus fluctuations determine effective 
2d behavior. The symmetry of isotropic vector TDKM 
is SU (2) and corresponding algebra is non-abclian. 

One may mention several more physical systems, where 
the scalar and/or planar TDQM is useful. One of such 
models is the big quantum dot with charge fluctua- 
tions accompanied by longitudinal and transversal spin 
fluctuations^. The class of Gaussian ensembles corre- 
sponds to infinite-range correlations in the charge and 
spin sectors of the model. Both spin and charge inter- 
actions contain stochastic component^, leaving a room 
for original formulation of the Keldysh model. The 
gauge field theory, based on functional bosonization be- 
ing formulated in the time-domain, opens a possibil- 
ity of stochastic treatment of dynamic processes. As is 
shown ir>i^, the transverse spin correlation function for 
anisotropic spin exchange contains both short-time and 
long-time correlation parts. While the short-time (white 
noise) correlations dominate away from the Stoner in- 
stability, the (infinitely) long-time correlations become 
important as one approaches the regime of strong fluctu- 
ations of the magnetization. The long-time part of the 
model is equivalent to the planar TDKM. 

Another interesting object is the optical superlattice 
consisting of biased double wells^. The bias is random, 
but the number of atoms in the same in all TLS in this ex- 
perimental setup due to the " interaction blockade" . One 
may expect that the well population in these TLS could 
be stochastized in accordance with Fig. 2J provided the 
Keldysh-type fluctuations ([5j or ([6]) with long relaxation 
times were realized experimentally. 

We thank Boris Altshuler, Jan von Delft, Yuri 
Galperin, Yuval Gefen, Jean Richcrt and Wilhclm Zw- 
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